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The magnetostatic interaction between two oppositely charged transverse domain walls (DWs) 
in adjacent Permalloy nanowires is experimentally demonstrated. The dependence of the pinning 
strength on wire separation is investigated for distances between 13 and 125 nm, and depinning 
fields up to 93 Oe are measured. The results can be described fully by considering the interaction 
between the full complex distribution of magnetic charge within rigid, isolated DWs. This suggests 
the DW internal structure is not appreciably disturbed by the pinning potential, and that they 
remain rigid although the pinning strength is significant. This work demonstrates the possibility of 
non-contact DW trapping without DW perturbation and full continuous fiexibility of the pinning 
potential type and strength. The consequence of the interaction on DW based data storage schemes 
is evaluated. 

PACS numbers: 75.75.+a, 75.60.Ch, 85.70.Kh, 75.60.Jk 
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Quantifying the interaction between magnetic domain 
walls (DWs) in closely packed networks of ferromagnetic 
nanowires is of vital importance to recently proposed DW 
based logic and data storage schemes 0,0, as these inter- 
actions could potentially compromise coherent DW prop- 
agation and correct device function. In addition to these 
schemes, DWs in nanowires have also been suggested for 
use in a wide range of applications, such as atom trap- 
ping for quantum information processing Q. Further- 
more, the fundamental properties of DWs themselves are 
attracting great interest, and DWs can be considered not 
only as transition regions between oppositely magnetized 
domains, but as quasi-particles Q. Their equilibrium 
dynamic properties [1, 0, H, [1] and interactions with ar- 
tificial defects [l^l are being intensively studied. The 
ability to hold DWs at defined positions is required for a 
wide range of spin torque experiments, where current is 
used to depin the DWs; a continuously variable pinning 
strength which does not alter the current path is highly 
desirable in such cases [llj . The possibility of a well de- 
fined, tunable pinning potential where the DW acts as a 
truly rigid quasi-particle is also appealing, for example, 
in current induced resonance experiments [l^ . [isj . The 
effect of the presence of a transverse DW (TDW) in a 
nearby structure on the magnetic state of a ferromag- 
netic ring has been reported fl3, HBl) and quantified in 
terms of an additional local field due to the TDW stray 
field, but no quantitative analysis of the full pinning po- 
tential created by a DW on another DW has been made. 

In this paper, we experimentally study the interaction 
between two oppositely charged DWs travelling in two 
parallel Permalloy (Py) nanowires and its dependence on 
wire separation. The separations probed in this investi- 
gation are below the dimensions of the DW (~ lOOnm) 
itself and so are within the near-field limit. By consid- 
ering the full, rigid magnetostatic charge distribution of 
an isolated DW we are able to reproduce fully the exper- 
imental results. This suggests the internal structure of 
the DW is not appreciably perturbed by the interaction 



for the geometries studied. This work opens up the pos- 
sibility of continuously tailorable, remote DW pinning, 
which does not affect the DW internal structure. 




FIG. 1: a) SEM image of a typical structure showing two 
lOfim long 'U' shaped Py nanowires separated by an edge 
to edge distance d — 100 nm. The asterisks indicate the 
positions where the MOKE laser spot is focused, b) Schematic 
of oppositely charged TDWs created in 'A' corners moving, 
under field H, along wires separated by d. 



100 nm wide, 8 nm thick, thermally deposited Py 
nanowires were fabricated using electron beam lithog- 
raphy (EBL) and a lift-off method. Fig. [T^ shows a 
scanning electron microscope (SEM) image of a typical 
structure: two U-shaped nanowires with an edge to edge 
separation, d, between 125 nm and 13 nm, and vertical 
arm displacement, w, equal to 0.5, 1 or 1.5 /im. Magnetic 
field sequences were applied in the plane of the sample 
at a frequency of 1 Hz using a quadrupole electromag- 
net. Two experiments were performed. The first was to 
create two oppositely charged DWs in the nearest pair of 
corners {Ai and A2). In this case (see Fig. [TJj), the DWs 
must pass within close proximity of one another before 
switching the remaining portion of each wire. The at- 
tractive interaction of the oppositely charged DWs will 
therefore directly influence the field required to separate 
them and switch the rest of the wire. In the second case, 
DWs are created in corners Bi and B2 and so are able to 
switch the central portion of each wire before meeting. 
This second setup acts as a control experiment to observe 
switching without any DW-DW interaction. 

To create DWs a.t A (B) corners a 45° field pulse was 
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applied in the +x, +y (—a;, +y) direction, giving an up- 
wards tail to tail (TT) DW at corner Ai (Bi) and upward 
head to head (HH) DW at corner A2 (-02)- The field was 
then decreased to zero at 45°, and a negative (positive) 
X field was applied to move the two DWs towards one 
another. The sequence was repeated with all fields re- 
versed to create a downward HH DW at Ai (Bi) and a 
downward TT DW aX A2 (S2), moving them towards one 
another with a positive (negative) x field. 

High sensitivity, spatially resolved magneto-optical 
Kerr effect (MOKE) measurements were used to probe 
the magnetization of the central portion of each wire 
(marked with *). Inset Fig. [5] is an example of the ob- 
served hysteresis loop for the switching of a nanowire. 
The 60 Oe switching field is due to the pinning from the 
DW-DW interaction. The switching fields of portions *i 2 
of the nanowires were measured as a function of d (mea- 
sured using SEM images). The results for 90 structures 
are shown in Fig. [2] Switching of the nanowires in the 
control experiment (DWs created in corners B) occurs 
at an average propagation field Hp of 19 ± 4 Oe (grey 
line). Fields H^i and iJ£)2 at which *i and *2 reverse 
in the case where the DWs can interact (DWs created at 
corners A) are seen to increase as d decreases. The mean 
difference between Hm and H02 is only 0.1 Oe with a 
standard deviation of 2.1 Oe; therefore only the average 
value Hu for each structure has been plotted (open cir- 
cles). We see a monotonic decrease of Hd as d increases 
from a maximum depinning field of 93 Oe measured in 
wires with d = 13 nm down to ~ 190e (~ Hp) for d '~ 90 
nm. Above d ~ 90 nm, the dependence changes as Ho re- 
mains approximately constant at Hp. At these values of 
d the dominant pinning mechanism is due to the inherent 
roughness of the nanowires and so the interaction of the 
DWs is no longer measured. No observable dependence 
was found on w. Separate studies, not presented here, 
ruled out a localised chan ge in wire geometry causing 
the observed DW pinning [l^, suggesting the observed 
Hp) is solely due to the interaction between the two DWs. 
This is also supported by the fact Hpi and Hp)2 are near 
identical for a pair of wires. The narrow distribution of 
Hp (as measured from control experiment) shows we are 
not measuring a statistical variation of Hp along these 
wires. 

To account for the DW-DW interaction dependence on 
d, we describe it in terms of the interaction between mag- 
netic charges, where the magnetic charge density pM is 
defined as Pm — —po^-M and M is the local mag- 
netization. We expect a model which considers the full 
DW charge distribution within the wire to match exper- 
imentally determined data if the two DWs remain rigid 
throughout the interaction and we account for finite tem- 
perature. OOMMF simulations [17] (3.5 x 3.5 x 4nm^ cell 
size, Ms = 800 kA/m, A = 13 x 10"^^ J/m and damping 
constant a = 0.5) were used to calculate the magnetiza- 
tion configuration of the transverse DW (TDW) which is 
stable for our wire geometry (100 nm wide, 8 nm thick) 
[1^. From this the charge distribution of the full 
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FIG. 2: Open circles: experimental dependence of DW trans- 
mission field Hd on wire separation d. Grey line, average wire 
propagation field Hp. Modelled depinning fields with a fit- 
ting scale factor /3 = 0.54, considering either a full rigid DW 
charge distribution (solid line), line-line interaction (dashed) 
or point-point interaction (dotted). Inset: Example MOKE 
observed hysteresis loop of a nanowire. Switching occurs at 
Hd, in this case 60 Oe, due to depinning from the DW-DW 
potential well. 



DW [20| was obtained. The charge density is very inho- 
mogeneously distributed, reflecting the triangular shape 
of TDWs [1^, with a positive surface charge along the 
wide edge of the DW, and a smaller negative charge along 
the narrow edge. 

The interaction energy, U, between two extended 
charge distributions, Di and D2, is given by: 
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(1) 



where is the distance between the two charges qi and 
Qj. In the case where the two distributions are associ- 
ated with oppositely charged rigid DWs travelling in two 
adjacent wires separated by d (see Fig. [Ub), [/ is a po- 
tential well, and depends on their lateral separation x. 
The change in Zeeman energy, SUz caused by separating 
the DWs by 5x under a fleld H is 5Uz = -~2MsHS6x, 
where Ms is the saturation magnetization and S is the 
cross sectional area of the wire. The potential landscape 
U{x) is sheared by a gradient equal to —2MsHS, and the 
field required for DW depinning is therefore given by: 
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Ma 



The solid line in Fig. [5] displays (depinning fields 

for the full charge distribution) as a function of d, fitted 
to the experimental data (by least mean squares) using a 
fitting scale factor (3 = 0.54. Fitting to the data is excel- 
lent up to d ~ 90 nm after which, as mentioned above. 
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the inherent wire roughness becomes the dominant pin- 
ning mechanism. 

This modelhng does not include the effect of temper- 
ature; however, DW depinning is a thermaUy assisted 
process. At room temperature, depinning fields have pre- 
viously been found to be between 40 and 65% of the zero 
temperature simulated values [2ll|. Ref. (2^ experimen- 
tally shows similar reductions between measurements at 
4.2 K and 300 K. The reduction in our modelling is well 
within the limits of these values. Our experimental re- 
sults are therefore in very good agreement with a full 
rigid charge model, suggesting the interacting DWs can 
be considered as rigid quasi-particles with identical in- 
ternal magnetic structure to that of a free, isolated DW. 

The far-field approximation of the charge distribu- 
tion to a point-like volume charge is the simplest an- 
alytical model possible and may agree with the near- 
field experimental data depending on the homogene- 
ity of the charge distribution. In order to test this, 
we consider a monopole containing the full DW charge 
Qbw = ^fiQMgS. A charge is located at the centre of 
each wire along the axis of symmetry of the DW, w/2 
away from the wide edge where w is the wire width. In 
this case, H^{d) is given by: 
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The resulting vs d curve is shown Fig. [2] (lower dotted 
line), using (3 = 0.54. The model is seen to underestimate 
Hjj over the region of d investigated due to the artificial 
concentration of the charge at a single point, showing we 
cannot make a far-field approximation in this region of 
d values and highlighting the inhomogeneity of the DW 
charge distribution. 
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FIG. 3; Potential well obtained when laterally separating a 
pair of oppositely charged DWs by x for the case d = 30 
nm. Solid line: obtained considering the DW as a finite, rigid 
distribution of magnetostatic charges. Circles: obtained using 
OOMMF simulations. 

Without using a detailed knowledge of the charge dis- 
tribution we could also attempt to form an analytical 
model based upon the geometry of the DW. TDWs have 
a triangular shape [18] as in Fig. [T]3 with a large amount 
of the magnetization along the wide edge aligned in the 



y direction. The charge due to the discontinuity of this 
magnetization across the wire edge will be significant and 
as the widest sides of the DWs are nearest we would ex- 
pect this to be the overwhelming contribution to the in- 
teraction. In order to test this, we consider a uniform 
line of magnetic surface charge, magnitude a per unit 
length and length 2L, with a x 2L = Qdw is used. The 
interaction energy between two such line charges of 
opposite sign, separated by (x, d) is found by solving 
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-L \/{xi - X2)^ + d'^ 



■ dxi 



■ dx-2 



(4) 

Hj^ for is shown FigO (upper dashed line) for the case 
where 2L is taken as the DW width and (3 — 0.54. For 
this the DW width used is the FWHM of the y compo- 
nent of magnetization along the wire centre in micromag- 
netic simulations (85 nm for our wire geometry). This 
attempt to use a more realistic analytical model with- 
out considering the detailed DW charge distribution is 
seen to overestimate the interaction strength. A more 
complete analysis of the distribution in our wire geom- 
etry shows only 48% of the DW charge is located along 
the wide edge, and has a Gaussian shape. Clearly then, 
positioning the full Qdw along the wire edge will over- 
estimate the interaction strength. We conclude therefore 
that only by considering the near-field full charge dis- 
tribution of the DW is it possible to closely match the 
experimentally observed magnetostatic interaction. 

Micromagnetic simulations were performed (parame- 
ters as above), modelling the situation as in Fig. [T]b. The 
resulting depinning field dependence on d is not shown, 
but is in good agreement with both experimental and 
near-field charge modelling (see Fig. [3]) and provides fur- 
ther evidence that little significant distortion occurs to 
the DW shape. 
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FIG. 4: Example of stray field based pinning creating an 
asymmetric potential well, a) Geometrical setup, b) Result- 
ing asymmetric potential well experienced by the DW. 

Using one DW to pin another is a specific example 
of pinning by stray fields with little distortion to DW 
shape. It should therefore be possible to create a pinning 
potential of arbitrary shape with careful use of localised 
stray fields. Fig. |4] shows an example of using multi- 
ple magnetic wire ends at different positions to create an 
asymmetric pinning potential. The depinning fields at 
K from each side of the well are 52 and 168 Oe respec- 
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tively, which will act as a DW diode or ratchet similar to 
the geometrical protrusions in Refs. [23,[2^. By modify- 
ing the geometrical layout of such wires pinning strength 
and well shape can be varied. Using such a tailorable 
potential that does not alter wire or DW shape in spin 
torque and current induced resonant experiments will al- 
low much clearer interpretation of current based effects. 

The interaction between DWs in neighboring 
nanowires is of great technological importance for 
DW based logic and data storage devices. Full charge 
modelling of DWs in nanowires of more technologically 
relevant dimensions (45 nm wide, 8nm thick, 45 nm 
spacing), shows DW depinning at 71 Oe (OOMMF 
simulations predict depinning at 67.5 ±2.5 Oe) at K. 
This large interaction could result in cross talk errors 
between adjacent data channels. 

In conclusion, we have experimentally studied the in- 
teraction between two DWs of opposite charge moving in 
parallel magnetic nanowires. The attractive interaction 



is found to result in depinning fields as high as 93 Oe, 
measured in wires with separations of 13 nm. Modelling 
using the full rigid magnetic charge distributions of iso- 
lated DWs reproduces the experimental results extremely 
well, suggesting that the DWs remain unperturbed by the 
interaction. Simple test analytical models that do not use 
the full charge distribution did not reproduce the exper- 
imental results, showing that consideration of the inho- 
mogeneity of the charge is essential to understand near- 
field DW interactions. Pinning fields at K of 71 Oe 
are expected in the technologically relevant dimensions 
of DW based data storage and logic devices. This study 
experimentally demonstrates remote pinning using mag- 
netostatic interaction with minimum perturbation to the 
DWs, and the possibility of continuous tailoring of pin- 
ning strength with no alteration to the nanowire. This 
work was supported by the Corrigan Scholarship and the 
European Community under the Sixth Framework Pro- 
gramme SPINSWITCH (MRTN-CT-2006-035327). 
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a) Micromagnetic simulation of transverse 
domain wall in lOOnm wide, 8nm thick 
permalloy nanowire (cell size 3.5x3.5x4nm^ 
M3 = 800kA/m, A = 1 3x1 0"''J/m). 

b) Magnetostatic charge distribution obtained 
from diverging magnetization in a). 

c) Magnetic surface charge per unit length 
along widest edge of domain wall for specific 
geometry shown. 



